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The one-dimensional harmonic oscillator in a box problem is possibly the simplest example of 
a two-mode system. This system has two exactly solvable limits, the harmonic oscillator and a 
particle in a (one-dimensional) box. Each of the two limits has a characteristic spectral structure 
describing the two different excitation modes of the system. Near each of these limits, one can 
use perturbation theory to achieve an accurate description of the eigenstates. Away from the 
exact limits, however, one has to carry out a matrix diagonalization because the basis-state 
mixing that occurs is typically too large to be reproduced in any other way. An alternative to 
casting the problem in terms of one or the other basis set consists of using an “oblique” basis 
that uses both sets. Through a study of this alternative in this one-dimensional problem, we are 
able to illustrate practical solutions and infer the applicability of the concept for more complex 
systems, such as in the study of complex nuclei where oblique-basis calculations have been successful. 
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I. INTRODUCTION 

The understanding of a physical system is closely 
linked to how well one can determine its eigenstates. 
Typically a set of basis states that works well in one limit, 
fails in another. And more general methods, such as vari¬ 
ational schemes, perturbation theory, or fixed-basis ma¬ 
trix diagonalizations typically begin with a reasonable 
Hamiltonian and some appropriate set of basis states that 
yield a good description of the system. 

When applying perturbation theory, one is usually con¬ 
cerned with a small perturbation of an exactly solvable 
Hamiltonian system. However, there are many exam¬ 
ples when the Hamiltonian has more than a single ex¬ 
actly solvable limit. This is a common situation when 
a dynamical symmetry group is used in the construc¬ 
tion of the Hamiltonian. A simple example is the hy¬ 
drogen atom in an external magnetic field: With increas¬ 
ing field strength, particularly for magnetic fields exceed¬ 
ing a so-called critical value of 2.35 x 10® T, the system 
changes from the spherical symmetry of the Goulomb 
problem to the cylindrical symmetry of the diamagnetic 
Hamiltonianii Another example that occurs in a variety 
of condensed matter problems is a particle confined to 
two dimensions in an external magnetic fields In nuclear 
physics, the Interacting Boson Model classifies many nu¬ 
clei according to one of three dynamical symmetries^ 
But, what should be done if the system is nowhere near 
any of the exact limits? In these situations, the prob¬ 
lem may be approached better by using states associ¬ 
ated with all the appropriate nearby limits. This set of 
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states will form an “oblique” (mixed-mode) basis for the 
calculation!^ In general, such a basis is non-orthogonal 
and may even be over complete. Nevertheless, as recent 
studies demonstrate, such oblique bases have merit. In 
this paper, we use a pedagogically simple problem to il¬ 
lustrate the oblique basis approach. 

The one-dimensional harmonic oscillator in a one¬ 
dimensional box has been used to illustrate different 
aspects of mixing generated by two interactions. Bar¬ 
ton, Bray, and Mackane used the model to study the ef¬ 
fects of distant boundaries on the energy levels of a one¬ 
dimensional quantum system^ Studies have also been 
done for the cylindrically symmetric system of a three- 
dimensional harmonic oscillator between two impenetra¬ 
ble wallsHowever, these studies did not discuss the 
bi-modal nature of the problem. Some authors have gen¬ 
eralized the one-dimensional harmonic oscillator prob¬ 
lem by introducing time-dependent parameters in the 
Hamiltonian^ to pass between the two limiting cases 
of a free particle and the harmonic oscillator solutions. 
The infinite square well and the harmonic oscillator have 
been considered as two limiting cases of a power-law po¬ 
tential within the context of wave packet collapses and 
revivals Recent research in modified uncertainty re¬ 
lations has also shown related and interesting behavior lii 

In this paper, we demonstrate the concept of an oblique 
basis approach by considering the simple two-mode sys¬ 
tem of the one-dimensional harmonic oscillator in a box. 
First, in Sec. QH we discuss the concept as well as the ex¬ 
actly solvable limits of this toy model. A qualitative dis¬ 
cussion of the expected spectrum of the one-dimensional 
harmonic oscillator in a box is given in Sec. cni along 
with an example spectrum and quantitative estimates. 
In Sec. nYi some specific problems related to the struc¬ 
ture of the Hilbert space are addressed. Sec. m contains 
specific toy model calculations and results as well as a 
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discussion of a new interesting behavior, similar to that 
observed in nuclear structure studiesi^ The discussion 
in Sec.0is focused on a quasi-perturbative behavior and 
a coherent structure within the “strong mixing” regioni^ 
where the system is far from any exact limit. Our conclu¬ 
sions are given in Sec. EDand suggested student exercises 
in Sec. IVIII 


II. HARMONIC OSCILLATOR IN A 
ONE-DIMENSIONAL BOX 

Let us start with an abstract two-mode system. For 
simplicity, we assume that the Hamiltonian of the system 
has two exactly solvable limits: 

H={l-X)Ho + XHi. (1) 

Clearly this is set up so that H Hq in the limit A ^ 0 
and H ^ Hi when A ^ 1. In the vicinity of these two 
limits, one can use standard perturbation theory for one 
Hamiltonian perturbed by the other^ Usually, some¬ 
where in between these two limits there is a critical value 
of A that is related to the strongly mixed regime of the 
system. This value of A could be anywhere in the inter¬ 
val (0,1), a convenient choice being A ~ i. Sometimes, 
a further symmetry breaking Hamiltonian H 2 can be ex¬ 
plicitly introduced by adding A(1 — X)H 2 to H. 

In Eq. (E), the variable A has been introduced to sim¬ 
plify the discussion. In general, there will be more than 
just one such parameter in the Hamiltonian^ Often the 
exactly solvable limits are described as hypersurfaces in 
the full parameter space. It could even be that there 
are three or more exactly solvable limits. For example, 
the Interacting Boson Model (IBM)^ for nuclear spectra 
has three exactly solvable limits>i^ii& Another example 
with three exactly solvable limits is the commonly used 
schematic interaction with non-degenerate single-particle 
energies (si), pairing {P^P) two-body interactions^ 
and quadrupole-quadrupoleiS (Q ■ Q) two-body interac¬ 
tions: 

H = e,N, - GP+P -xQ-Q- 

Here, we consider what is perhaps the simplest two¬ 
mode system that shares the essential features of such 
problems while remaining pedagogically instructive. The 
Hamiltonian of a one-dimensional harmonic oscillator in 
a one-dimensional bo»iS of size 2L has the form: 

= + + ( 2 ) 

where VL{q) is the confining potential taking the value 
zero for \q\ < L and 00 for jgl > L, and w is the oscillator 
frequency. This system has two exactly solvable limits. 

The first limit of the toy model in Eq. © is w = 0 
when it reduces to a free particle in a one-dimensional 
box of size 2L, 

Ho = + VL{q). (3) 


The eigenvectors and energies are labeled by n = 0,1,... 
and given by 


^n{q) 


Er, 



((n-b l)f ■§) if n is even 
((n -b l)f I;) if n is odd 


( 4 ) 



This limit corresponds to extreme nuclear matter when 
the short range nuclear force can be described as an effec¬ 
tive interaction represented by a square-well potentiali^S 
We can think of this limit as the one-dimensional equiv¬ 
alent of a three-dimensional model where nucleons are 
confined within a finite volume of space representing the 
nucleus. Recently such effective potentials for the Bohr 
Hamiltonian have been used to introduce symmetries in 
the critical point of quantum phase transitions^ 

The other exactly solvable limit of the toy model in 
Eq. 121), when L —> 00 , is the harmonic oscillator in one 
dimension, 


Hi = - 

2m 


-q ■ 


In dimensionless coordinates. 


( 5 ) 


q 



pV mhuj, 


we have 


Hi = {p^ -b q^) . 

The eigenvectors and energies are labeled by u = 0,1,... 
and are given by 


En = /zw(n-bi), b=J—, 

\ 2 J V mto 

where Hn are the Hermite polynomials. This limit is 
essentially the harmonic oscillator model for nuclei. 

In a one-dimensional toy model, the anharmonic oscil¬ 
lator with a quartic anharmonicity would be the appro¬ 
priate counterpart of the Sp{6, R) shell modelSS since the 
quadrupole-quadrupole interaction Q ■ Q goes as ~ r'^ and 
Q connects harmonic oscillator shells with like parity. If 
we restrict the model space to only one harmonic oscilla¬ 
tor shell, then we can use the algebraic quadrupole mo¬ 
ment Q of Elliott^Si because within a single shell Q is the 
same as Thus, for single shell studies, it is appropri¬ 
ate to consider the one-dimensional harmonic oscillator 
as representative of the SU (3) shell model for nuclei. 
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III. SPECTRAL STRUCTURE AT DIFFERENT 
ENERGY SCALES 

Often in physics the spectrum of a system has dif¬ 
ferent characteristics over different energy regimes. This 
usually reflects the existence of different excitation modes 
of the system. For the toy model Hamiltonian in Eq. 
we can define three spectral regions: 

• Spectrum of a particle in a one-dimensional box as 
in Eq. with quadratic dependence on n {E^ ~ 

• Spectrum of the one-dimensional harmonic oscil¬ 
lator as in Eq. @ with linear dependence on n 

• Intermediate spectrum that is neither of the above 
two types. 


This analysis is confirmed by the numerical calculations 
shown in Fig. 13 where the perturbed particle in a box 
spectrum provides a good description for n > 6 for the 
case oi m = h — 2 L/'k — 1 and w = 4. Actually, the 
agreement extends to lower n as is often the case, pertur¬ 
bation theory seemingly yielding valid results well past 
its expected region of validity. Note that the first or¬ 
der corrections are already close to the limiting constant 
value of On the other hand, first-order pertur¬ 

bation clearly fails for the ground state. Indeed, an ear¬ 
lier studjiiS of the inadequacy of fixed basis calculations 
showed that adequate convergence for the ground state 
when uj is large requires a large number of basis states 
of the one-dimensional box. This is related to the fact 
that for large values of w, or equivalently large L values, 
a large number of particle-in-a-box wave functions (sin 
and cos) are needed to obtain the correct behavior (ex¬ 
ponential fall off) of the low-energy harmonic oscillator 
wave functions in the classically forbidden zone. 



Particle in Harmonic Particle in a ID box subject to 

ID box Oscillator harmonic oscillator potential. 



FIG. 1: Two-mode toy system consisting of a particle in a 
one-dimensional box subject to a central harmonic oscillator 
restoring force. 

As shown in Fig.Q] one expects to see the particle in 
a box spectrum at high energies. These energies corre¬ 
spond to the box boundaries dominating over the har¬ 
monic oscillator potential. In this regime, one can use 
standard perturbation theory to calculate the energy for 
a particle in a box perturbed by a harmonic oscillator po¬ 
tential. It can be shown that perturbation theory gives 
better results for higher energy levels. For n —> oo, the 
first correction [5El^) approaches a constant value: 

SEl = (l - (5^;^) - 

Using UO+i - EO » (n |U| n), one estimates that pertur¬ 
bation calculations are valid when 


FIG. 2: Exact energies of a two-mode system with m = h = 
2L/7r = 1 and w = 4 compared to the spectrum of the one¬ 
dimensional harmonic oscillator (left), spectrum of the free 
particle in a ID box (right), and spectrum as calculated within 
a first order perturbation theory of a free particle in a ID box 
perturbed by a ID HO potential. The lowest three eigenener- 
gies of the two-mode system nearly coincide with the ID HO 
eigenenergies, while higher energy states are better described 
as perturbations of the other limit of a free particle in a ID 
box. 

The intermediate spectrum should be observed when 
the harmonic oscillator turning points coincide with the 
walls of the box. Therefore, the critical energy that sep¬ 
arates the two extreme spectral structures is given by 


Notice that the constant energy shift tow^L^/ 6 in the 
energy of the high energy levels SE^^^ is one-third of 
the critical energy (Ec/3). 


n > 


3fi^7r2 
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At low energies, where the classical turning points of 
the oscillator lie far from the boundaries, we expect the 
spectrum to coincide essentially with that of the oscillator 
as shown in Fig.|21 The number of such nearly harmonic 
oscillator states is easily estimated using 


Ec > E. 


HO 


^HO 


1 mujL^ 

2 h 


1 

2' 


(9) 


There is also a compatible number of levels, usually larger 
than below the Ec corresponding to a free particle 

in a box. 


Ec > El^ ^ nlR^ = 


2 mujL^ 


- 1 . 


( 10 ) 


However, these states are mixed by the harmonic oscilla¬ 
tor potential. 

Using the ratio of the ground state energies, 
Eq^/Eq^ = AmojL '^together with Eq. @ and 
Eq. lUni) . the following spectral situations apply: 


• Eor (f)^ < , there are levels below Ec corre¬ 

sponding to the harmonic oscillator and the parti¬ 
cle in a box such that Eq'^ > Eq^. The oscillator 
levels dominate the low energy spectrum. 

• Eor ^ < (f)^, there are only the ground 

states Eq^ and Eq^ below Ec and Eq^ > Eq^. 

• Eor 1 < there is only the ground state 

of the harmonic oscillator E^^ below Ec- 

• For 0 < — < 1, perturbation theory of a parti¬ 

cle in a box should be applicable for all levels. 

The dimensionless parameter jS = thus plays a role 

similar to the parameter A in Eq. O- In our abstract case 
of Eq. (^3), the two limits of H are at A = 0 and A = 1 
with strong mixing at A = 1/2. In the case above, the 
two limits are /? = 0 and /? = oo with strong mixing when 
1 < P < (7r/2)^. In this respect, we can formally relate 
P to X using an expression of the form A = P/{P + Pc) 
where Pc is the value of P in the strong mixing region. 
For example, one may chose Pc = 7r/2. 

Consider as a numerical illustration of the two-mode 
spectra the case ot m = h = 1, L = 7r/2 and w = 4 shown 
in Fig. [3 With these parameters, Eq. gives = 
4.53. Thus one should see no more than 4 equidistant 
states. Indeed, in Fig.|3 there are four equidistant energy 
levels that correspond to a harmonic oscillator spectrum. 

With respect to the critical energy Ec, there is a more 
explicit classification of the spectral structure: 


• Perturbed particle in a one-dimensional box spec¬ 
trum for energies E ^ Ec such that Eq. 0 holds, 

• One-dimensional harmonic oscillator spectrum in 
Eq. ® for energies Ec^ E such that Eq. @ holds, 

• Intermediate spectrum for energies E k, Ec- 


IV. HILBERT SPACE OF THE BASIS WAVE 
FUNCTIONS 

Before discussing the toy model using an oblique ba¬ 
sis, it is instructive to discuss briefly the harmonic oscil¬ 
lator problem in Eq. 0 using the wave functions for a 
free particle in a one-dimensional box as in Eq. 0 ; and 
vice versa, solving the problem of a free particle in a one¬ 
dimensional box in Eq. 0 using the wave functions for 
a particle in the harmonic oscillator potential given in 
Eq. 0- 

Due to the different domains of the wave functions, 
there are some specific problems that need to be ad¬ 
dressed. For example, using wave functions for a free 
particle in a one-dimensional box to solve the pure har¬ 
monic oscillator problem may not be appropriate espe¬ 
cially for high energy states E ^ Ec- This is because 
all such functions vanish outside the box (see Fig. 0) un¬ 
like the oscillator wave functions, especially as the energy 
increases. The converse is also problematic because os¬ 
cillator functions that are non-zero outside the box will 
lead to infinite energy. 



FIG. 3: Spreading of the wave functions: harmonic oscillator 
wave functions spread outside the harmonic oscillator poten¬ 
tial into the classically forbidden region; particle in a box wave 
functions are zero at and outside of the box boundary. 

The influence of the boundary conditions on the prop¬ 
erties of a quantum-mechanical system has been rec¬ 
ognized from the dawn of quantum mechanics. It is 
well known that some problems with seemingly sepa¬ 
rable Hamiltonians may re-couple due to the boundary 
conditions^ Some recent studies on the problem of con¬ 
fined one-dimensional systems using equations for rele¬ 
vant cut-off functions have been discussed by Barton, 
Bray, and Mackanei^ Their method has been further de¬ 
veloped in a more general setting by Berman^ Other au¬ 
thors aim at variational procedures using simple cut-off 
functions2Si2i or derive asymptotic estimates for multi¬ 
particle systems using the Kirkwood-Buckingham varia- 
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tional methodi^ Somewhat different approaches focus on 
shape-invariant potentials and use supersymmetric part¬ 
ner potentials to derive energy shifts and approximate 
wave functionSfSS as well as dependence of the ground- 
state energy on sample size^ In the next few paragraphs, 
we discuss the structure of the relevant Hilbert spaces 
when confinement is present. 


A. Harmonic Oscillator in the One-Dimensional 
Box Basis 

Consider the harmonic oscillator in Eq. using the 
wave functions for a free particle in a one-dimensional box 
in Eq. 0 . There are no practical difhculties for energies 
E Ec as defined in Eq. 0 as long as the turning 
points of the oscillator are sufficiently deep into the box 
(sufficiently far from the walls of the box). However, for 
energies E ^ Ec, the basis wave functions are localized 
only on the interval [—L,L] and thus cannot provide the 
necessary spread over the width of the potential (Fig. 0. 
This situation would be appropriate for the toy model in 
Eq. 0 but not for the pure harmonic oscillator problem 
in Eq. 0. 

One simple solution to the spreading problem is to con¬ 
tinue the basis wave functions using periodicity. This way 
the necessary spread of the basis wave functions can be 
achieved and the new basis will stay orthogonal but must 
be re-normalized. If one continues the wave functions to 
infinity, normalization will require Dirac delta functions 
but for continuation on a finite interval, the functions 
can be normalized to unity as usual. However, these ba¬ 
sis wave functions do not decay to zero in the classically 
forbidden zone. This means that a significant number of 
basis wave functions will be needed to account for the ap¬ 
propriate behavior within the classically forbidden zone. 

Another alternative is to change the support domain 
corresponding to non-zero values of the function by 
stretching or squeezing, accomplished through a scaling 
of the argument of the basis wave functions, x xan/L. 
This way the support becomes [—L,L\ —> [—a„,Q!„]. 
Here, is a scale factor for the n-th basis wave function 
in Eq. 0 , estimated either from the width of the har¬ 
monic oscillator potential, or determined by variational 
minimization. Either way, the new set of basis functions 
will be non-orthogonal. In general, there may even be 
a linear dependence. However, for the basis functions 
discussed here, linear dependence may not appear due 
to the different number of nodes for each wave function. 
The number of nodes (zeros) is not changed under the 
re-scaling procedure. While the potential width scal¬ 
ing is simpler, its applicability is more limited than the 
variationally-determined one which can be extended to 
more general situations^ 


B. Particle in a Box in the Harmonic-Oscillator 
Basis 

When the choice of the basis is not made with appro¬ 
priate care, an operator that should be Hermitian may 
become non-hermitian. In nuclear physics, although this 
is unlikely for a finite shell-model calculation using an oc¬ 
cupation number representationj^i it is an obstacle when 
one wishes to use a hard core potential and a harmonic- 
oscillator basis 

Suppose we want to solve the problem of a free particle 
in a one-dimensional box [—A, L] as given in Eq. 0 using 
the harmonic oscillator wave functions in Eq. 0 . The 
first thing to do is to change the inner product of the 
wave functions: 

/ oo pL 

f*{x)g{x)dx^ / f*{x)g{x)dx. (11) 

-OO J —L 

Then, it is immediately clear that the set of previously 
orthonormal harmonic oscillator wave functions 
in Eq. 0 will lose their orthonormality and even their 
linear independence. The set of functions (q) with 
support domain restricted to [—L,L] and denoted by 
'f'rt(<Z; [—L, L]) become linearly dependent if L is so small 
that there is more than one function dt„(g; [—L, L]) with 
the same number of nodes within [—L,L]. While this 
linear dependence can be handled, the real problem is 
the loss of hermiticity of the physically significant op¬ 
erators. Neither the variational nor the potential-width 
wave function scaling will help to cure the loss of her¬ 
miticity of operators, such as those of linear momentum 
and energy. This non-hermiticity is due to the behav¬ 
ior of the basis states at the boundary, mainly the non¬ 
vanishing of the wave functions at —L and L. 

To understand the loss of hermiticity, we look at the 
off-diagonal matrix elements of the momentum operator 
{p = After some trivial manipulations, we have: 

(«'rn,P^'„) = [p'Hm. «'n) - ih (^-^ ((Z)^'„ (?)) | • 

It is clear from the above expression that hermiticity will 
be maintained only when all of the basis functions are 
zero at the boundary of the interval [-L, L\. Actually, the 
necessary condition is that the wave functions have the 
same value at ±L; they need to be zero only for the case 
of an infinite potential wall. 

In our simple example, all operators are built from the 
momentum operator p and position operator q. Thus, 
in order to ensure hermiticity, it is sufficient to make 
sure that p is Hermitian which requires the basis wave 
functions to vanish at the boundaries —L and L. For 
this purpose, one can look at the nodes of each basis 
wave function and scale it so that its outer nodes are at 
the boundary points. From the nodal structure of the 
harmonic oscillator wave functions, it is clear that the 
first two wave functions (d>o and d>i) cannot be used 
since they have fewer than two nodes. Since the physical 
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requirement that the wave functions have to be zero at 
the boundary is the cornerstone in quantizing the free 
particle in a one-dimensional box as in Eq. 0 , it is not 
surprising that the nodally adjusted harmonic oscillator 
wave functions are very close to the exact wave functions 
for the free particle in a one-dimensional box as shown 
in Fig. ^ 


(a) (b) (c) 


FIG. 4: Harmonic-oscillator trial wave functions (dark gray) 
adjusted with respect to the one-dimensional box prob¬ 
lem: (a) adjusted according to the potential width = 

^ uin = ^ {1 + 2n), (b) nodally adjusted (first three 
are deliberately phase shifted), (c) boundary adjusted using 
®(g) ^ ®(g) - (L) (1 -b q/L)/2 - ^ (-L) (1 - q/L)/2. The 

exact wave functions (light gray) for a particle in a box are 
zero at ±1, as clearly seen in the (a) graphs. 

In general, calculating the nodes of a function may 
become very complicated. To avoid problems with find¬ 
ing the roots, one can evaluate the wave function at the 
boundary points, then shift the wave function by a con¬ 
stant to get zeros at the boundary, 'I'(( 7 ) ^ fE'(( 7 ) — (L). 

This idea works well for even parity wave functions, but 
has to be generalized for odd parity by adding a linear 
term, f['(g) ^ ^'(g) — (Ik (L) /L)q. Thus, for a general 
function, we can have: 'l'(g) —> 'i>{q) — {l+q/L)'^ (L) /2 — 
(1 — q/L)'i> {—L) /2. In Fig.^lwe have shown some of the 
resulting wave functions. Notice that this procedure gives 
a new wave function df that has good behavior inside the 
interval [—T, L] and grows linearly with q outside the in¬ 
terval \—L,L\. This contrasts with the behavior of the 
cut-off function f{q) obtained by Barton et al^ where the 
function /(g) has an L/q singularity at the origin (g = 0). 
The use of a cut-off function to enforce boundary condi¬ 
tions has been studied by Barton et a& and BermaiiS^ 
and provides an interesting integral equation for the cut¬ 
off function. On the other hand, a simple cut-off function 
supplemented by a variational method seems to be very 
effect ive 

An alternative, more involved construction can be ex¬ 
plored which relies on the Lanczos algorithm. This algo¬ 
rithm is an iterative process that uses the Hamiltonian 


of the system to generate a new basis state from the pre¬ 
vious state at each iteration. By using the boundary 
matching process above, one can set up and successfully 
run a modification of the usual Lanczos algorithm^S to 
solve for the few lowest eigenvectors of the free particle 
in a one-dimensional box through an arbitrary, but rea¬ 
sonable, choice of the initial wave function^ The major 
modification is to project every new function, dtn+i = 
H'i’n, into the appropriate Hilbert space and subtract 
the components along any previous basis vectors. Only 
then should one attempt to evaluate the matrix elements 
of H related to the new basis vector that is clearly within 
the correct Hilbert space. This way, one has to dou¬ 
ble the number of scalar product operations compared 
to the usual Lanczos algorithm where the matrix ele¬ 
ments of H are calculated along with the complete re- 
orthogonalization of the basis vectors. 


C. Oblique Basis for the Two-Mode System 

In the previous two subsections we saw that a basis 
that is appropriate for one mode of our model system, 
harmonic oscillator in a box, is not appropriate for the 
other mode. If we desire a description in the critical 
mixed region characterized by Ec as shown in Fig. ^ 
a combination of the two basis sets seems appropriate. 
This is referred to as an oblique basis, stemming from a 
geometrical analogy. A two-dimensional space is usually 
described in terms of the Cartesian (x, y) axes or, indeed, 
in terms of any other orthogonal pair obtained by rotat¬ 
ing those axes. While such orthogonal choices are more 
convenient and familiar, any two axes, as long as they are 
not linearly dependent, also suffice to describe the space. 
Such a choice constitutes an “oblique” pair. Similarly, 
when we mix both harmonic oscillator and particle in a 
box states, we have an oblique basis. 

In using such an oblique basis, there are two main prob¬ 
lems to be addressed in order to have a proper Hilbert 
space of our quantum mechanical system. First, we have 
to make sure that any set of states that are derived from 
harmonic oscillator functions satisfies the boundary con¬ 
ditions of the problem. For the particle in a box states the 
boundary conditions are satisfied by construction. We 
discussed already a few possible ways to construct states 
with the correct boundary conditions. An interesting ad¬ 
ditional method, suggested by one of our referees, would 
be to use Hermite functions with non-integer index. For 
such functions with index between n and n -I- I, the po¬ 
sition of the outer node is correspondingly between the 
outer nodes of the n-th and (n-l-l)-th Hermite functions. 
Second, after the chosen set of functions has been mod¬ 
ified to satisfy the boundary conditions, the orthonor¬ 
mality of these functions would most likely be destroyed. 
Even if the two basis sets are orthonormal by themselves, 
they would not be orthonormal as a whole and may even 
be linearly dependent. While this might seem a serious 
technical problem, it has a well known solution through 
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re-orthonormalization of the basis or by proceeding with 
a generalized eigenvalue problem4i 

Our oblique basis consists of modified harmonic oscil¬ 
lator (MHO) basis states that satisfy the boundary con¬ 
ditions along with basis states of a free particle in a box 
(BOX). To obtain the wave functions and eigenenergies 
we solve the generalized eigenvalue problem within this 
oblique basis. Schematically, these basis vectors and their 
overlap matrix can be represented in the following way: 


basis vectors 
overlap matrix 
hamiltonian 


E = 
0 = 
H = 


Be, : box — basis \ 
Ci : mho — basis J 

i b Ei t - Cq 


\El+ fi J Hij = 




HryS Hfy 


ajS \ _ 

-^i/3 Elij J 

I EIi)Qxxbox El^Qxxmho 
\ EljjihQxbox EljYiboxmho 


where a = 1,..., dim(box - basis) and * = 1,..., dim(mho 
- basis). 

In these notation, the eigenvalue problem: 


Hv = Ev 


with 


V = u“ea + V^Bi 

takes the form: 



which is a generalized eigenvalue problem. 


V. DISCUSSION OF THE TOY MODEL 
CALCULATIONS AND RESULTS 

Despite the simplicity of the toy model in Eq. @, 
the harmonic oscillator in a box exhibits some of the es¬ 
sential characteristics of a more complex system. Some 
of our interest lies in problems associated with the use of 
fixed-basis calculations. In particular, one such prob¬ 
lem is the slow convergence of the calculations If 
one can implement an exact arithmetic, one may not 
worry too much about the slow convergence when enough 
time, storage, and other computer resources are provided. 
However, numerical calculations are plagued with numer¬ 
ical errors so that a calculation that converges slowly 
may be compromised by accumulated numerical error. 
Of course, having the correct space of functions is of the 
utmost importance in any calculation. 

Having considered the main problems one may face 
in studying the simple toy model in Eq. , we con¬ 
tinue our discussion with the spectrum for the case of 
m = h = 2 L/tt = 1 and w = 4. As one can see in 


Fig.m the first three energy levels are equidistant and 
almost coincide with the harmonic oscillator levels as ex¬ 
pected from Eq. ®. For these states, the wave func¬ 
tions are also essentially the harmonic oscillator wave 
functions. The intermediate spectrum is almost miss¬ 
ing. Above Ec, the spectrum is that of a free particle in 
a one-dimensional box perturbed by the harmonic oscil¬ 
lator potential. We find that an oblique-basis calculation 
reproduces the lowest eight energy levels using 14 ba¬ 
sis functions, seven nodally adjusted harmonic oscillator 
states and seven states of a free particle in a box. In 
contrast, a fixed-basis calculation, using only the wave 
functions of a free particle in a one-dimensional box, re¬ 
quires 18 basis states. 

Due to the simplicity of the toy model, it does not 
appear as if the oblique-basis calculation has a big nu¬ 
merical advantage over calculations using the fixed basis 
of the box wave functions. There are two main reasons 
for this: (1) there is a critical energy that separates 
the two modes, (2) the spectrum above Ec has a nice 
regular structure. 

The regular structure above Ec results in a very favor¬ 
able situation for the usual fixed-basis calculations since 
the dimension of the space needed to obtain the n-th 
eigenvalue grows as n + a. The parameter a is relatively 
small and does not change much in a particular region 
of interest. For example, the w = 16 calculations need 
only a = 15 extra basis vectors when calculating any 
of the eigenvectors up to the hundredth. The relatively 
constant value of a can be understood by considering the 
harmonic oscillator potential as an interaction that cre¬ 
ates excitations out of the n-th unperturbed box state. 
Therefore, a is the number of box states with energies in 
the interval and -|- uj‘^/2 (<i)„| |$„) where is 

the n-th unperturbed box state energy. There is a rapid 
de-coupling of the higher energy states from any finite 
excitation process that starts out of the n-th state. This 
de-coupling is due to the increasing energy spacing of the 
box spectrum which results in a limited number of states 
mixed by the harmonic oscillator potential. Using the 
upper limit Ec/i on SE^, one can estimate a ~ 

The sharp separation of the two modes also allows for 
a safe use of the harmonic oscillator states without any 
rescaling. This is especially true when to is very large 
since then the low energy states are naturally localized 
within the box. Therefore, instead of diagonalizing the 
Hamiltonian in a box basis, one can just use the harmonic 
oscillator wave functions. 

Fig. 0 shows the absolute deviation {AE = — 

j^esumate'^ of the Calculated energy spectrum for the case 
of w = 16, L = 7r/2, h = m = l. Here, stands 

for one of the three energy estimates one can make: the 
harmonic oscillator E^^, particle in a one-dimensional 
box E^^, and the first order perturbation theory esti¬ 
mate considering the harmonic oscillator potential as a 
perturbation, {E^^+cu'^/2 ($„| |4>„)). There are about 
19 states that match the harmonic oscillator spectrum 
which is consistent with the expected value from Eq. 0. 
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FIG. 5: Absolute deviations of variously calculated energies 
from the exact energy eigenvalues for w = 16, L = 7r/2, h = 
m = 1 as a function of n. Circles represent deviation of 
the exact energy eigenvalue from the corresponding harmonic 
oscillator eigenvalue (Ai? = — E^^), the diamonds 

are the corresponding deviation from the energy spectrum of 
a particle in a box (Ai? = — El^), and the squares are 

the first-order perturbation theory estimates. 


After the n = 20 level, perturbation theory gives increas¬ 
ingly better results for the energy eigenvalues. Fig. El 
shows the relative deviation (1 — 



FIG. 6: AsinFig.Elbut for relative deviations from the exact 
energy eigenvalues of the three calculations. 

Note that perturbation theory is valid, as expected, for 
high energy states determined by Eq. 0 . For the high 
energy spectrum, the harmonic oscillator potential acts 
as a small perturbation. Thus the first-order corrections 
in the energy and the wave function are small. Fig. [ 7 | 
shows that the main component of the 105*^ exact wave 
function comes from the 105*^ box wave function, as it 
should in a region of small perturbations. 

For low energy states, perturbation theory around 


FIG. 7: Non-zero components of the 105th exact eigenvector 
in the basis of a free particle in a one-dimensional box. Pa¬ 
rameters of the Hamiltonian are uj = 1%, L = 7r/2, h = m — 1. 


the box states is not appropriate, the harmonic oscil¬ 
lator states being closer to the true states in this region. 
Specifically, for m = h = 2L/7r = 1 and to = 16, the first 
ten states are essentially the harmonic oscillator states 
to a very high accuracy. The next ten states have still 
high overlaps with the corresponding harmonic oscilla¬ 
tor wave functions. For example, starting from 0.999999 
at the tenth state, the overlaps go down to 0.880755 at 
the twentieth state. After that the overlaps get small 
very quickly. Fig. El shows the structure of the third ex¬ 
act eigenvector when expanded in the box basis. Notice 
that the third box wave function is almost missing. An 
explanation lies in the structure of the harmonic oscilla¬ 
tor functions, which are essentially exact in this region. 
Upon projecting these functions in Eq. onto the box 
functions in Eq. 0 , the results can be obtained in closed 
analytical form. The integrand consists of three factors, 
an even power of q, a Gaussian and a cosine. While the 
oscillations of the Hermite polynomial are of varying am¬ 
plitude, the cosine has evenly spaced nodes and antinodes 
of equal amplitude. As a result, cancellations can take 
place between successive terms in the integrand. For the 
third oscillator function (second of even parity), its one 
pair of nodes is reflected in the dip seen at n = 3 in Fig. El 
Higher oscillator functions with more pairs of nodes can 
display correspondingly more such dips in the projected 
squared amplitudes. 

This pattern of having a small projection of the exact 
wave function along the corresponding box wave function 
continues to persist into the transition region. This may 
seem unexpected considering the fact that the first order 
estimates of the energy levels are relatively good. Notice 
that in our two examples of w = 4 (Fig. 0 and w = 16 
(Fig. 0, the first order corrections in the transition re¬ 
gion are already close to the limiting constant value of 
even though the corresponding box wave func- 
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FIG. 8: Non-zero components of the third harmonic oscillator 
(even parity) eigenvector as expanded in the basis of a free 
particle in a one-dimensional box. Parameters of the Hamil¬ 
tonian are uj = 16, L = vr/2, ft = m = 1. 


tions are not present at all in the exact wave function as 
shown in Fig.lHl 

From the results in these graphs, it seems that the 
transition region is absent since first-order perturbation 
theory becomes valid immediately after the breakdown of 
the harmonic oscillator spectrum. That first-order per¬ 
turbation theory gives good estimates for the energy lev¬ 
els in this transition region is a manifestation of coherent 
behavior. What actually happens in this region is a co¬ 
herent mixing of box states by the harmonic oscillator 
potential in the sense of a quasi-symmetryi^4^2i24 This 
coherent mixing is illustrated in Fig. where one can 
see that the histograms for a few consecutive states are 
very similar. In this sense we say that the corresponding 
particle in a box states are coherently mixed. A more 
precise and detailed discussion of coherent structure, be¬ 
havior, and quasi-structures can be found elsewherei^ 


VI. CONCLUSIONS 

In summary, we have studied the simplest two-mode 
system consisting of a one dimensional harmonic oscil¬ 
lator in a one-dimensional box. Depending on the pa¬ 
rameters of the two exactly solvable limits, one observes 
various spectral structures. There is a clear coherent mix¬ 
ing in the transition region. It is remarkable that such 
a simple system can exhibit coherent behavior similar 
to the one observed in nuclei. There is clear advantage 
to using an oblique-basis set which includes both oscil¬ 
lator and particle in a box states. This allows one to 
use the correct wave functions in the relevant low and 
high energy regimes relative to E^. Taking into account 
the importance of the relevant energy scale of a prob¬ 
lem and the wave function localization with respect to 
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FIG. 9: Goherent structure with respect to the non-zero com¬ 
ponents of the 25th, 27th and 29th exact eigenvector in the 
basis of a free particle in a one-dimensional box. Parameters 
of the Hamiltonian are lj = 16, L = 'k/2, ft = m = 1. 


the range of the potential, the oblique-basis method can 
be extended beyond the idea of using two or more or¬ 
thonormal basis sets. Specifically, one can consider a 
variationally-improved basis set by starting with some 
initially guessed basis states. In the occupation number 
representation (Fock space), which is often used in the 
nuclear shell model^ this variationally-improved basis 
method seems inapplicable. But the method is of gen¬ 
eral interest because of its possible relevance to multi¬ 
shell ab-initio nuclear physics, atomic physics, and gen¬ 
eral quantum mechanical calculations. The method may 
also be related to renormalization-type techniquesi^ 


VII. STUDENT EXERCISES 

1. Section II views the two Hamiltonians in Eq. © 
and Eq. © as limits of Eq. for suitable choices 
of parameters. Express this feature instead in the 
form of Eq. © with suitable definition and choice 
of A and appropriate modification of u in Eq. ©• 

2. Treating the harmonic oscillator potential as a per- 
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turbation, work out the first order correction to the 
energies of particle in a box states. Hence, verify 

Eq. 0. 

3. Project the third harmonic oscillator wave function 
(second even parity state with two nodes) onto the 
wave functions of a particle in a box to verify the 
structure shown in Fig. |H1 Integrals involved have 
products of powers of q, a Gaussian and a cosine in 
q, and you may take integration limits to ±oo. 

4. As a variational counterpart of the oblique basis 
concept in this paper, the product of the ground 
state functions of a particle in a box and an oscil¬ 
lator may be chosen as a trial wave function. With 
such a form, calculate through the Rayleigh-Ritz 
variational principle the ground state energy of the 


system. 
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